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Precise measurements of the production and decay offer a window for searches for new physics. In this pa- 
properties of the heaviest known fermion, the top quark, per we measure the top-antitop quark pair (tt) produc- 
provide important tests of the standard model (SM) and tion cross section and compare it with the SM prediction, 
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extract the top quark pole mass from this measurement 
and search for new physics in top quark decays analyzing 
ratios of the tt cross sections measured in different decay 
channels. 

The inclusive tt production cross section (<r t j) is mea- 
sured in different tt decay channels assuming SM branch- 
ing fractions. The comparison of the results to pre- 
dictions in next-to-leading order perturbative quantum 
chromodynamics (QCD), including higher order soft 
gluon resummations [J, 0, S 0] , yields a direct test of 
the SM. Ratios of c t f measured in different final states 
are particularly sensitive to non-SM particles that may 
appear in top quark decays, especially if the boson in 
the decay is not a SM W boson. An example is the 
decay into a charged Higgs boson (t — > H + b), which, 
as predicted in some models [5], can compete with the 
SM decay t — ► W + b. Additionally, many experimental 
uncertainties cancel in the ratios. Furthermore, since a t i 
depends on the mass of the top quark (m t ), it can be used 
to extract mt- Such measurement is less accurate than 
direct mass measurements, but provides complementary 
information with different experimental and theoretical 
uncertainties. 

Within the SM, each quark of the tt pair is expected 
to decay nearly 100% of the times into a W boson and 
a b quark @. W bosons can decay hadronically into qq' 
pairs or leptonically into ev ei [iv^ and tv t with the r 
in turn decaying onto an electron, a muon, or hadrons, 
and associated neutrinos. If one of the W bosons decays 
hadronically while the other one produces a direct elec- 
tron or muon or a secondary electron or muon from r 
decay, the final state is referred to as the £+jets (or £]) 
channel. If both W bosons decay leptonically, this leads 
to a dilepton final state containing a pair of electrons, a 
pair of muons, or an electron and a muon (the 11 chan- 
nel), or a hadronically decaying tau accompanied either 
by an electron or a muon (the t£ channel). 

Measurements of the individual tt cross sections in ££ 
and t£ channels using about 1 fb _1 of pp data from the 
DO detector at the Fermilab Tevatron collider at ^fs = 
1.96 TeV are available in Ref. Q- In the £+jets channel, 
we use the same selection and background estimation as 
in Ref. Q, but a slightly larger dataset and a unified 
treatment of systematic uncertainties with the ££ and t£ 
channels. We provide a brief summary of the event se- 
lection and analysis procedures below. 

In each final state we select data samples enriched in 
tt events by requiring one or two isolated high transverse 
momentum (pr) leptons for the ^+jets or ££ channel re- 
spectively. At least two high pt jets are required for ££ 
and t£ events, and at least three for £+jets events. Fur- 
ther, in all but the e/i channel, large transverse missing 
energy {$t) is required to account for the large trans- 
verse momenta of neutrinos from W boson or r lepton 
decays. In the e/x final state, a requirement on the sum 
of the pr of the highest pr (leading) lepton and the two 
leading jets is imposed instead. In the /i/j, channel, the 
$t requirement is supplemented with a requirement on 



the significance of the $t measurement, estimated from 
the pr of muons and jets, and their expected resolutions. 
Additional criteria are applied on the invariant mass of 
the two opposite charge leptons of the same flavor in the 
ee and fj.fi channels to reduce the dominant background 
from Z/-f* — > £ + £~ events. In the £+jets and t£ chan- 
nels we require a minimum azimuthal angle separation 
between the J^t vector and the lepton px, A(J)(£,^t), to 
reduce background from multijet events, where jets are 
misidentified as electron, muon or r. Details of lepton, 
jet and $t identification are provided in Refs. @, Hoj - 
The final selection in these channels demands at least one 
identified b jet via a neural- network based algorithm [Til ] . 
In the £+jets channel we separate events with one or > 2 
b-tagged jets due to their different signal over background 
ratio and systematic uncertainties. 

To simplify the combination and extraction of cross 
section ratios, all channels are constructed to be mutu- 
ally exclusive. In particular, events with two identified 
leptons are excluded from the £+jets selection, and all t£ 
candidates are removed from the rest of the channels. 

The compositions of the samples in the ^+jets, ££ and 
t£ channels are shown in Table HI I-F+jets production 
dominates the background for the £+jets events, while 
multijet production is the most important background 
in the t£ channel. Background in the ££ channels comes 
mainly from Z+jets production. In the ££ channel, con- 
tributions from IF+jets production are part of the mul- 
tijet background. The smaller contribution from dibo- 
son production is included in the category labeled "other 
background" . This category also includes the contri- 
bution from single top quark production in the £+jets 
and t£ channels. The signal, VF+jets and Z+jets back- 
grounds are simulated using ALPGEN [l^] for the matrix 
element calculation and PYTHIA [l3| for parton showering 
and hadronization. Diboson and single top backgrounds 
are simulated with the PYTHIA and SINGLETOP [lj] gen- 
erators, respectively. We estimate the multijet back- 
ground from the control data samples. The difference 
in the ratio of tt and W+jets events in the e+jets and 
/j+jets final states is the result of the larger efficiency and 
misidentified lepton rate in the e+jets channel compen- 
sating for the lower lepton acceptance (|7y| < 1.1) com- 
pared to the /z+jets channel (|?7| < 2.0). In addition, the 
wider rapidity distribution of the M^+jets events com- 
pared to tt, events increases the W+jets background con- 
tribution in the yu.+jets channel. 

To calculate the combined cross section, we define a 
joint likelihood function as the product of Poisson proba- 
bilities for the 14 disjoint subsamples, as listed in TableQ] 
Fourteen additional Poisson terms constrain the multijet 
background in the £+jcts and t£ channels. In particular, 
for the re and t/x channels, the multijet background is 
determined by counting events with an electron or muon 
and associated t of the same electric charge, introducing 
a corresponding Poisson term per channel. In the £+jets 
channel, we estimate the multijet background separately 
for each of the eight subchannels by using corresponding 
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TABLE I: Expected numbers of background and signal events for a t i — 8.18 pb, observed numbers of data events and measured 
<r t t at top mass of 170 GeV/c 2 . Quoted uncertainties include both statistical and systematic uncertainties, added in quadrature. 



Channel 


Luminosity(pb i ) 


W+jets 


Z+jets 


Multijet Other bkg 


tt 


Total 


Observed 


era (pb) 


e+jets (3 jets, 1 b tag) 


1038 


53A±$» 


e-o±i:S 




11.4+i'^ 


8i.7ts : ? 




183 


8.0611$ 


^t+jets (3 jets, 1 b tag) 


996 


59-2+56 




9.7^28 


9.5±i : i 


59.0±^ 


143.9 + ^i 


133 


6-43+J; 2 , 2 


e+jets (3 jets, > 2 b tags) 


1038 


5.0±g-| 




2.7+^ 

— U.o 


2-4+U1 

— U.4 


30.7^3.9 


4i.5±*:i 


40 


7.781^1 


^t+jets (3 jets, > 2 b tags) 


996 


5.8±g;i 


o.7±g:3 


l o +0 - 3 


9 1 +0.3 
** A — 0.3 


90 0+3.4 

Zo.o_3 2 


33 .5 g g 


31 




e+jets (> 4 jets, 1 b tag) 


1038 


8.51^7 


9 +U.b 


7.9+!-" 

■ — l.U 


3.0+^ 

— U . 5 


81.6^91 


103.3t 7 ;g 


113 


9.38+i;r 2 


/^+jets (> 4 jets, 1 b tag) 


996 


13.6±i;f 


2.5 q g 


0+ 00 


9 ,+0.4 


65.9±?;1 


84.3t 6 i 


99 


10.44l 2 ; 7 J 


e+jets (> 4 jets, > 2 b tags) 


1038 


1 o+ u - 3 

i- u -o.3 


2+ u + 


1 l+ u + 
L -- L -0.1 


9 +u ' 2 
u.a_ 2 


' -6.0 


44.9^-0 


30 


c 19+1. w 

d - lz -1.28 


/i+jets (> 4 jets, > 2 b tags) 


996 




3+ 01 


0+ 00 


0.7l° :i 


35.6+ 5 5 ? 


38.2+1* 


34 


7.60l 2 ; 7 J 


ee 


1074 




n o+U.5 
■ 4 -°-0.5 


6+ u ' 4 
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control data samples [15(. Four additional terms arise 
from applying this same method in evaluating the multi- 
jet background before b tagging. 

Each systematic uncertainty is included in the like- 
lihood function through one free "nuisance" parame- 
ter [l5|. Each of these parameters is represented by a 
Gaussian probability density function with zero mean 
and a standard deviation of one; all are allowed to float 
in the maximization of the likelihood function, thereby 
changing the central value of the measured a t i. Correla- 
tions are taken into account by using the same nuisance 
parameter for a common source of systematic uncertainty 
in different channels scaled by the corresponding stan- 
dard deviation each individual channel. Thus, the like- 
lihood function to be maximized is represented by the 
product 

14 14 K 

C = Y[P(n u mi ) x Y[V(n jt mj) x JJSD i)fc X Q{v k ; 0, 1) , 

i=l j = l k=l 

where P(n,m) is the Poisson probability to observe n 
events given the expectation of m events. The predicted 
number of events in each channel is the sum of the pre- 
dicted background and expected events, which depends 
on cr t j. In the product, i runs over the subsamples and 
j runs over the multijet background subsamples. The 
Gaussian distributions SD^-x^i/^; 0, 1) describe the sys- 
tematic uncertainties, K is the total number of indepen- 
dent sources of systematic uncertainty, are the indi- 
vidual nuisance parameters, and SDj/. is one standard 
deviation for the source of uncertainty k in subsample i. 

Systematic uncertainties on the measured a t i are eval- 
uated from sources that include electron and muon iden- 
tification; r and jet identification and energy calibra- 



tion; 6-jet identification; modeling of triggers, signal and 
background; and integrated luminosity. All these uncer- 
tainties are treated as fully correlated among channels 
and between signal and background. Systematic uncer- 
tainties arising from limited statistics of data or Monte 
Carlo samples used in estimating signal or backgrounds 
are considered to be uncorrelated. A detailed discussion 
on systematic uncertainties can be found in Refs. @>[8[. 
Table HI1 shows a breakdown of uncertainties on the com- 
bined cross section. We evaluate the effect from each 
source by setting all uncertainties to zero except the one 
in question and redoing the likelihood maximization with 
respect to only the corresponding nuisance parameter. 
Since the method allows each uncertainty to change the 
central value, the total uncertainty on o t i differs slightly 
from the quadratic sum of the statistical and individ- 
ual systematic uncertainties. The total systematic un- 
certainty on a t i exceeds the statistical contribution. The 
luminosity uncertainty of 6.1% which enters into the esti- 
mation of the majority of the backgrounds and the lumi- 
nosity measurement of the selected samples is the domi- 
nant source of systematic uncertainty. 

Table Hill summarizes the individual a t i measurements 
for the individual channels, as well as some of their com- 
binations. Within uncertainties, all measurements are 
consistent with each other. The combined cross section 
for £+jets, It and rt final states for a top quark mass of 
170 GeV/c 2 is evaluated to be 

<r tf = 8.18l8;|? pb, (2) 

in agreement with theoretical predictions [l], [1, H, H| . The 
uncertainty is comparable to the one on the cross section 
combination from different methods in the £+jets channel 
performed by DO 8]. The observed number of events 
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FIG. 1: (a) Expected and observed numbers of events versus channel, used in measuring the combined a t i. The dashed band 
around the prediction indicates the total uncertainty. Upper limits on B(t — ^ H + b) for (b) tauonic and (c) leptophobic H + 
decays. The yellow band shows the ±1 standard deviation band around the expected limit. 



TABLE II: Summary of uncertainties on the combined <j t t- 



Source 




(Pb) 


Statistical 


+0.47 


-0.46 


Lepton identification 


+0.15 


-0.14 


Tau identification 


+0.02 


-0.02 


Jet identification 


+0.11 


-0.11 


Jet energy scale 


+0.19 


-0.16 


Tau energy scale 


+0.02 


-0.02 


Trigger modeling 


+0.11 


-0.07 


b jet identification 


+0.34 


-0.32 


Signal modeling 


+0.17 


-0.15 


Background estimation 


+0.14 


-0.14 


Multijet background 


+0.12 


-0.12 


Luminosity 


+0.56 


-0.48 


Other 


+0.15 


-0.14 



Total systematic uncertainty +0.78 —0.69 



TABLE III: Summary of measured a t i in different channels 
for m t = 170 GeV/c 2 . 



Channel 


<?ti (pb) 


^+jets 
11 [7] 
^+jets and It 

H w 

l+jets, U and rl 


8.46± ;° 9 7 

10+0.99 

o-lo_0.87 
7 77 +2.90 

1 ■ II -2.4,7 
1O+0.98 
O-J-O-0.87 



in the different channels is compared to the sum of the 
background and combined tt signal in Fig. QJa). 

We compute ratios R a of measured cross sections, 
I&M = <tII/ct% and R T J /U ^ = ag/vg*", by gener- 
ating pseudo-datasets in the numerator and denomina- 
tor in order to take into account the correlation between 
systematic uncertaintes. -c| lannel represent the measured 
cross sections in the corresponding channel. The pseudo- 
datasets are created by varying the number of signal and 
background events around the expected number accord- 
ing to Poisson probabilities. All independent sources of 



systematic uncertainties are varied within a Gaussian dis- 
tribution. Although the individual channels considered 
are exclusive, each channel can receive signal contribu- 
tions from different tt decay modes. We calculate the con- 
tribution from dilepton events to the ^+jets final state as 
well as the contribution from dilepton and £+jets events 
to the t£ final states using the corresponding observed 
cross sections in the individual channels when generating 
pseudo-datasets. For each pseudo-dataset, we perform 
the maximization of Eq. [1] separately in the numerator 
and denominator, and divide the results. The central 
value is obtained from the mode of the distribution of 
R„, and the uncertainties are derived from the interval 
containing 68% of the pseudo-experiments. From these 
pseudo-experiments we obtain R f ^^ = 0.861q ^ and 
R ri/U-ej = o g7 +o.32 ^ which is cons i s tent with the SM 
expectation of R a = 1. 

Extensions of the SM, based on supersymmetry or 
grand unification Q, require the existence of additional 
Higgs multiplets beyond the Higgs doublet of the SM. 
Some of these models, such as the Two Higgs-Doublet 
Model or the Minimal Supersymmetric Standard Model, 
foresee the existence of physical degrees of freedom which 
can be associated with a charged scalar particle, the 
charged Higgs boson. If this charged Higgs boson is 
lighter than the top quark, it will appear in the top quark 
decays. We use the ratios to extract upper limits on 
the branching ratio B = B(t — * H + b). In particular, 
a charged Higgs boson decaying into a tau and a neu- 
trino (B(H + — * tv) = 1) results in more events in the tI 
channel, while fewer events appear in the li and £+jets 
final states compared to the SM prediction. In case of 
a leptophobic (B(H + — > cs) — 1) model, the number 
of dilepton events decreases faster than the number of 
£+jets events for increasing B(t — ► H + b). We therefore 
use R 1 ^^ to set limits on the leptophobic model, while 
jgt/u-ij . g use( ^ to searc h f or decays in which the charged 
Higgs bosons are assumed to decay exclusively to taus. 

To extract the limits, we generate pseudo-datasets as- 
suming different branching fractions B(t — > H + b). The 
signal for a charged Higgs boson is simulated using the 
PYTHIA Monte Carlo event generator [l3j], and includes 
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decays of ti — > W + bH~b and its charge conjugate (WH) 
and — > H + bH~b (HH). For a given branching frac- 
tion B, we calculate the expected number of ti events per 
final state, 

N t - t = [{l-B) 2 -e ww +2B{l-B)-e W H+B 2 -e HH ]<j tl L , (3) 

where e are the selection efficiencies for the different de- 
cays (WW refers to ti — > W + bW~b) and L is the in- 
tegrated luminosity. We add N t i to the expected back- 
ground and treat the sum as a new number of expected 
events in each channel. We then perform the likelihood 
maximization to extract a t i from these pseudo-data as 
if they contained only SM ti production. This provides 
distributions for the ratios of cross sections for each gen- 
erated B, which are compared to the observed ratio. We 
set limits on B by using the frequentist approach of Feld- 
man and Cousins [l6| . 

The observed and expected (i.e., for R a = 1) limits 
for the tauonic and the leptophobic charged Higgs boson 
models are shown in Figs.QJb) andQJc), respectively. In 
the tauonic model the upper 95% CL limits on B range 
from 15% to 40% for 80 GeV/c 2 < M H+ < 155 GeV/c 2 , 
improving the limits given in (l7| . For the leptophobic 
charged Higgs boson model, which is investigated here for 
the first time, the upper limit on B range between 48% 
and 57% for the same mass range. Although indirect 
bounds as those from the measured rate of b — ► S7 [HI 
appear stronger than the results from the direct search 
presented here, they can be invalidated by the presence 
of new physics contributions. 

The interpretation of the direct measurement of the 
top quark mass [6|, has become a subject of intense dis- 
cussion in terms of its renormalization scheme [19J . The 
extraction of this parameter from the measured cross sec- 
tion provides complementary information, with different 
sensitivity to theoretical and experimental uncertainties, 
relative to direct methods that rely on kinematic details 
of the top quark reconstruction. Simulated samples of ti 
events generated at different values of the top quark mass 
are used to estimate the signal acceptance. The resulting 
measurements of a t i are fitted as a function of m t 0]: 

Vtt( m t) = —[a+b(mt-m )+c(m t -mo) 2 +d(m t -mo) 3 ] 
m t 

(4) 

where cr t i and m t are in pb and GeV/c 2 , respectively, and 
rn = 170 GeV/c 2 [2(j. The dependence on the top mass 
is due to the mass dependence of the selection efficiencies. 

We compare this parameterization to a prediction in 
pure next-to- leading-order (NLO) QCD to a calcula- 
tion including NLO QCD and all higher-order soft-gluon 
resummations in next- to- leading logarithms (NLL) 0, 
to an approximation to the next-to-next-to-leading-order 
(NNLO) QCD cross section that includes all next-to- 
next-to-leading logarithms (NNLL) relevant in NNLO 
QCD [3j, and to a calculation that employs full kinemat- 
ics in the double differential cross section beyond NLL 
using the soft anomalous dimension matrix to calculate 
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FIG. 2: Experimental and theoretical [2,[2, HI °tt a s function 
of rat. The colored dashed lines represent the theoretical un- 
certainties due to the choice of the PDF and the renormaliza- 
tion and factorization scales. The point shows the measured 
combined a t i, the black dashed line the fit with Eq. [4] and the 
gray band the corresponding total experimental uncertainty. 

the soft-gluon contributions at NNLO J4j. Figure [5] shows 
the experimental and the theoretical [JH, Q ti cross sec- 
tions as a function of the top quark mass. 

Following the method of Refs. 0, @|, we extract the 
most probable top quark mass values and the 68% CL 
band. Since the theoretical predictions are performed in 
the pole mass scheme, this defines the extracted param- 
eter here. The results are given in Table llVl All values 
are in good agreement with the current world average of 
171.2 ±2.1 GeV/c 2 @. 

TABLE IV: Top quark mass with 68% CL region for different 
theoretical predictions of cr t j. Combined experimental and 
theoretical uncertainties are shown. 



Theoretical prediction 


mt (GeV/c") 


NLO [1] 


165.5t 6 5 i 


NLO+NLL [2] 


167.5; 


H5.8 

-5.6 


approximate NNLO [3] 


169.1; 


h5.9 
-5.2 


approximate NNLO [4] 


168.2; 


h5.9 

-5.4 



In summary, we have combined the tt cross section 
measurements in £+jets, tl and t£ channels to measure 
a ti = 8.18±g:lf pb for a top quark mass of 170 GeV/c 2 . 
For the first time, we have also calculated ratios of cross 
sections and interpreted them in terms of limits on non- 
standard model top quark decays into a charged Higgs 
boson. All results are in good agreement with the SM 
expectations. Finally, using different theoretical predic- 
tions given in the pole mass scheme, we have extracted 
the top quark mass from the combined a t i and have found 
the result to be consistent with the world average top 
quark mass [|| from direct measurements. 
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